Continuous time quantum walks in phase space 
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We formulate continuous time quantum walks (CTQW) in a discrete quantum mechanical phase space. We 
define and calculate the Wigner function (WF) and its marginal distributions for CTQWs on circles of arbitrary 
length N. The WF of the CTQW shows characteristic features in phase space. Revivals of the probability 
distributions found for continuous and for discrete quantum carpets do manifest themselves as characteristic 
patterns in phase space. 
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I. INTRODUCTION 



The study of quantum mechanical transport phenomena on 
discrete structures is of eminent interest. One particular as- 
pect is the crossover from quantum mechanical to classical 
behavior. Since classical physics is described in phase space 
and quantum mechanics in Hilbert space a unified picture 
is desired. This is provided, for instance, by the so-called 
Wiener function (WF) flU, which has remarkable proper- 
ties: It transforms the wave function of a quantum mechanical 
particle into a function living in a position-momentum space 
similar to the classical phase space. The WF is a real valued 
function and in this respect compares well with the classical 
probability density in phase space. However, it is not always 
positive. 

The concept of phase space functions has been widely used 
in Quantum Optics |3, 4] but also for describing electronic 
transport, see e.g. Here, the formal similarity to the 

classical Boltzmann distribution has been exploited. However, 
the phase spaces considered there are mostly continuous and 
infinite. This results in a simpler mathematical tractability of 
WFs than for discrete and for finite systems. Nevertheless, 
one can also define WFs for discrete systems 

Describing the transport by coined or continuous time 
quantum walks has been v ery successful over the last few 
years, for an overview see fol . The applicability of quan- 
tum walks reaches far beyond quantum computation; for in- 
stance in quantum optics also the (discrete) Talbot effect can 
be described by continuous time quantum walks 1 14, Q. An 
additional, but more abstract approach to transport processes 
is through the so-called quantum multibaker maps flU flr^l . 
Such maps were shown to exhibit (as a function of time) 
a quantum-classical crossover, where the crossover time is 
given by the inverse of Planck's constant. 

In the following we define the WF for a continuous time 
quantum walk on a one-dimensional discrete network of arbi- 
trary length N with periodic boundary conditions (PBC). We 
further show that the marginal distributions are correctly re- 
produced. Additionally, we define a long time average of the 
WF, which we compare to the classical limiting phase space 
distribution. Finally, we show how (partial) revivals of the 
probability distribution manifest themselves in phase space. 



II. CONTINUOUS TIME QUANTUM WALKS 

The quantum mechanical extension of a continuous time 
random walk (CTRW) on a network (graph) of connected 
nodes is called a continuous time quantum walk (CTQW). It 
is obtained by identifying the Hamiltonian of the system with 
the (classical) transfer matrix, H = T, see e.g. fTsL fl^l 
(we will set h = 1 in the following). The transfer matrix of 
the walk, T = (Ty), is related to the adjacency matrix A of 
the graph by T = — 7 A, where for simplicity we assume the 
transmission rate 7 of all bonds to be equal. The matrix A has 
as non-diagonal elements Ay the values —1 if nodes i and j 
of the graph are connected by a bond and otherwise. The di- 
agonal elements An of A equal the number of bonds /, which 
exit from node i. 

The basis vectors \j) associated with the nodes j span the 
whole accessible Hilbert space to be considered here. The 
time evolution of a state \j) starting at time t is given by 
\j\t) = V(t,t )\j), where U(t,t ) = exp[-zH(i - to)] is 
the quantum mechanical time evolution operator. 



III. WIGNER FUNCTIONS 
A. Definition 

The WF is a quasi-probability (in the sense that it can be- 
come negative) in the quantum mechanical phase space. For a 
d dimensional system, the (quantum mechanical) phase space 
is 2d dimensional. If the phase space is spanned by the con- 
tinuous variables X and K, the WF is given by Ijjy] 



W{X,K;t) = i /" dY . 



iKY 



{X-Y/2\p{t)\X + Y/2), 
(1) 

where p(t) is the density operator and thus for a pure state, 
p(t) = \i){t)){ij){t)\. Here, tp(X;t) = (X\ij)(t)) is the wave 
function of the particle. Integrating W(X, K; t) along lines in 
phase space gives marginal distributions, e.g., when integrat- 
ing along the if -axis one has, Illy], 



j dKW(X,K;t) 



MX;t)f 



(2) 



In the case of a discrete system, given, for instance, by N 
discrete positions on a network, which we choose to enumer- 
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ate as 0, 1, . . . , N — 1, the functions ip(x) are only defined for 
integer values of x = 0, 1, . . . , N — 1, and the form of Eq.Q 
has to be changed from an integral to a sum. There have 
been several attempts in doing so, see, for instance, iHH]. 
However, the definition of the discrete WFs might depend on 
whether the length N of the system is even or odd 110,201. 

In the following we proceed somewhat differently than pre- 
vious works and focus on discrete systems with PBC. For a 
one-dimensional system of length N (exemplified by a ring) 
this implies i/j(x) = ip(x ± rN) for all r € N. It follows that 
each and every one of the products ip*(x — y'; t)ip(x + y'\ t) 
is identical to (at least) one of the N forms ip* (x — y; t)ip(x + 
y; t), where y — 0, 1, . . . , N — 1. 

Now the WF has the form of a Fourier transform; a unique 
transformation of these N products requires N different k- 
values. These fc-values may evidently be chosen as k — 
2ttk/N, again having k = 0, 1, . . . , N— 1. We are thus led to 
propose for integer x and y the following discrete WF 



N-l 



w(x, M) = i£ eikv ^ x " y ' + y; *)• 



(3) 



y=o 



As a side remark we note that the rhs of Eq.Q stays in- 
variant when replacing N by mN (with m integer and m ^ 
0). Furthermore, also the choice of the summation interval 
0, 1, . . . , N — 1 is arbitrary; any N consecutive y-values will 
do. Given that the nature of the Hilbert space underlying 
Eq.0 differs from that of Eq.(|3}, there is no simple relation 
connecting the two equations. Here, one may recall the prob- 
lem of the normalization of plane waves in bounded (compact) 
and in infinite spaces. In the following we will use Eq.(|3} in 
order to compute the WFs. 



B. WFs and the Bloch ansatz 

In the following we consider only CTQWs with nearest 
neighbor steps. The Hamilton operator for such a CTQW on 
a finite one-dimensional network with PBC takes on a very 
simple form 



H|j) = 2\j) - \j - 1) - |j + 1) 



(4) 



where we have taken the transmission rate to be 7 = 1. The 
states \j), where j = 0, 1, . . . , N — 1, are the basis states as- 
sociated with the nodes j. Since here the CTQW describes 
a particle moving in a periodic potential we can employ the 
Bloch ansatz 1171 12111 . The time independent Schrodinger 
equation H|$g) = Eg\$>g) has eigenstates \&g) which are 
Bloch states. The PBC require that $g(N) = <I>e(0), where 
§e(x) — As before, this restricts the ^-values to 

9 = 2nn/N, where n = 0, 1, . . . , N - 1. The Bloch state 
can be expressed as a linear combination of the states 

|J>. 



I$6 
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N-l 



(5) 



In turn, the state | j) can be expressed by the Bloch states, i.e. 
by a Wannier function l2lll . 



3) = 



-iSj 



$6 



(6) 



where the sum runs over the N values for 9. From this it 
follows readily that the eigenvalues to Eq.@ are given by 
Eg = 2 — 2cos6>. The time evolved state then follows as 



\j;t) ^ e-* m |j) = 1 ^e-^e-^'l^), (7) 



Using the position basis means introducing ipj(x;t) = 
(x\j; t). Evidently, at t = one has ipj(x; 0) = (x\j) = S x j. 
Furthermore, with Eq.0 we have 

V>*(x';f)^(z;i) = (j;t\x')(x\j;t) 

6,6' 

1 e -H0-e') je i(,B e ,-E e )t e as e~ a ' s ' S X ,, V 5 X>1 



N 2 



LI' 



= J_ e - l(e - e ' ): > e l(Ee '~ Ee)t e iex e~ ie ' x ' . 

^ 6,6' 

Now, following Eq.(|3j, the WF is given by 



(8) 



W, 



(x,k;t) = 



x J2 e 

y 



i(6+6'+k)y 



(9) 



The sum over y can be written as NAg + g> + k where, as an 
extension of the usual Rronecker symbol, we introduce A„, 
such that A„ = 1 for n = (mod N) and A„ = else. 
Therefore, the WF for a CTQW on a one-dimensional network 
of N nodes with PBC reads 

1 N ^ 

Wj(x, n-t) - ^ ex P [ ~ i2 < 2n + K )( x ~ 3)/N] 



n=0 



x exp { - i2t[cos(2ir(n + n)/N) - cos(27rn/A^)]}, (10) 

where we have used that 9 = 2nn/N and k = 2itk/N. Fur- 
thermore, one may note from Eq.Q that 



ty/N 



j=0 



k k y 

xiP*{x-y;t)^(x + y;t) 
= H(x;t)\ 2 . (11) 

Note that the WF in Eq.( I10> does not change with time for 
k = 0, i.e. we have Wj(x,0;t) = M^(x,0;0) = Wj{x,0), 
wh ere Wj (x, n) is the long time average of the WF given in 
Sec llVCT 
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IV. WFS AND CTQWS 

In the following we consider WFs for CTQWs on circles 
of different length. For t — we start the CTQWs from one 
node. In particular, we take N — 100 (as an example of an 
even value) and N = 101 (as odd value); as starting site we 
choose j = 50. 



A. CTQWs in phase space 

In general, the WFs have a very complex structure. Fig- 
ures □ and display the WFs for CTQWs over cycles with 
N = 101 and N = 100, respectively, at t = 40. The peaks at 
(x, k) = (j, 0) for N = 101 and for N = 100 are clearly no- 
ticeable (black peaks in Figs[0and|2j; however, the remaining 
structure is quite complex. 




FIG. 1 : (Color online) 3D visualization of the WFs of a CTQW on a 
cycle of length N = 101 at t = 40. The initial node is at j = 50. 




FIG. 2: (Color online) Same as Figfl] for N = 100 at t = 40 with 
j = 50. 



Figure [5] shows a contour plot of the WF of a CTQW on a 
cycle of length N = 101 at different times. Note that at t = 
the WF is localized on the strip at the initial point j. At t = 1, 
the WF is still mostly localized about j. As time increases 
other sites get populated. On short time scales, the WF devel- 
ops a very regular structure in phase space, with "wavefronts" 
originating from the initial point j. Additionally, we also note 
"wavefronts" starting from the region opposite to the initial 
point, which are much weaker in amplitude. As time pro- 
gresses, these two types of waves start "interfering" with each 
other. 

In FigHwe plot the WF of a CTQW for N = 100. The 
structure of the WF is quite similar to the odd numbered net- 
work. Nonetheless, there are differences at larger times, visi- 
ble by comparing Figs|3je) and|3jf) to FigsOe) and|4jf). 




FIG. 3: (Color online) WFs of a CTQW on a cycle of length TV = 101 at times t= 1,10,20,40 [(a)-(d)] as well as t = 100,500 [(e),(f)]. The 
initial node is at j — 50. Red regions denote positive values of Wj(x, k; t), blue regions negative values and white regions values close to 0. 
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FIG. 4: (Color online) Same as Fig|3 for N = 100 and j = 50. 



However, although at long times the interference effects are 
quite intricate, typical patterns are still visible. At t = 500, 
we find less regularities in phase space for TV = 101 than 
for N — 100, reflecting the higher symmetry of CTQWs on 
even numbered networks. We note moreover, that the phase 
space pictures give at any time a much richer picture of the un- 
derlying dynamics than the transition probabilities \ijjj(x; t)\ 2 
alone. While the transition probabilities appear to be quite ir- 
regular 1 17], the WFs display regular patterns. As in the case 
of a particle in the box 12211 . also here we might view the en- 
suing structure, the "quantum carpets", to be woven by the 
discrete WFs. 

We finally note that in an infinite system the WFs have a 
much simpler structure, because there we do not face the prob- 
lem of distinct wave fronts running into opposite directions 
and interfering with each other. 



B. Marginal distributions 

In general, the marginal distribution J2 K Wj (x, k; t) is ob- 
tained from Eq.(I10> as 



1 N ~ 1 

Y,W 3 {x,K-t) = -J2ex.p[-i2nn(x-j)/N] 

K n— 

x exp[z2tcos(27rn/JV)] 

JV-l 

i 

X 



— exp[-i27r(/c + n)(x - j)/N] 

K=0 



x exp[-j2£cos(27r(ft + n)/N)]. (12) 

Since the system is periodic, the arguments (k + n) in the 
exponentials can be substituted by (k + n) = N — v, where 



v = 0,1,..., N - 1. This yields 

AT-l 



^2Wj{x,K;t) 



N 



(x—j)/N e —i2tcoa(2-jTv/N) 



(13) 

which is exactly what also follows, see Eg. ( II 1> . from calculat- 
ing \ip(x; t)\ 2 directly from the Bloch ansatz, see Eq.(|8}. The 
same relation was also highlighted in Ref.| 131 . 

Furthermore, we obtain the marginal distribution 
S K Wj' ( x ' K ' ^ or tne infinite system by taking N large (but 
fixed) in Ea. (ll3> . and setting 9 = 2-kv/N. We are then led to 



2tt 

— / d Q e iHx-j) -i2tco S ( 




(14) 



where J u (t) is the Bessel function of the first kind, see 
Eq.(9.1.21)ofRef.j2l. 

For the marginal distribution J2 x Wj(x, K;t) we also get 
from Eg. (Ti0t 



N-l 



)j/N 



x n=0 

x exp { - i2t[cos(27r(K + n)/N) - cos(27rn/iV)]} 
'l/N N odd and all k 

2/N N even and n even (15) 
N even and k odd, 

which is independent of t. Eq.(ll5> can be confirmed directly 
by taking the Fourier transform of \ip(x; t)\ 2 . The whole phase 
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space volume is normalized to one, as can be seen by sum- 
ming Eq.jl5> over all k, with k = 0, 1, . . . , N — 1. 



C. Long time averages 

Similar to the classical case, the CTQW will visit any point 
in phase space. However, due to the evolution with the unitary 
H, there is no definite limiting distribution. For the CTQW 
one may view the long time average of the transition proba- 
bility \ipj(x; t)\ 2 as a limiting distribution |24j|, i.e. 



= lim 



T 



(16) 



Employing this idea, we define the limiting WF as 



Wj(x, k) = lim 

T^oc 



- / dtWj{x,K',t). 



and for even N 



e _ I [2N- 2)/N 2 for x = j, j + N/2 
Xxj ~ 1 {N-2)/N 2 



else 



(21) 



which confirms our previous calculations 1 17]. 



D. Revivals in phase space 

CTQWs on large networks show only partial revivals fl7ll . 
These should occur at times t r }t N 2 /2tt. However, for small 
networks of sizes N = 1,2, 3, 4, 6 also complete revivals are 
possible, as was recently seen in waveguide arrays 11411 . The 
WFs are easily obtained from Ea. (llO> . The marginal distri- 
butions J2 K Wj(x, k; t) are, of course, correctly recovered. 
For instance, for N = 6, the first complete revival occurs for 
t = 2ir. This is readily seen, since for the WF at x = j we 
have Wj(x = j, k; t = 2n) = 1/6; thus the marginal distri- 
bution is Wj(x = j, K;t = 2ir) = 1, which also means, 



(17) using Ea.llTlll. that \tpj(x =j;t= 2tt)\ 2 = 1 



Upon integration of Eq. dlOi over t the only term which re- 
mains is the one for which cos[27t(k + n) /N] = cos[27rn/iV]. 

Some care is in order here, since for the long time average 
there are differences between even and odd N, also occurring 
in the limiting distribution Xxj IU7I1 . In the case of an odd N 
(superscript °) we get 



W°Ax,k) 




k / and any x 
k = and x = j 
else. 



(18) 



Note that for k ^ we require 2n = N — k. 

For an even N (superscript e ), the limiting WF for k = 
has contributions only for x — j and x = j + N/2. Each of 
the two contributions has weight 1/2. Thus the limiting WF 
reads 



1 /N 2 k/0 and any x 
W?(x,k) = { 1/2N k = and x = j, j + N/2 
else. 



(19) 



Note that also here we require that 2n = (N — k). 

The long time average of the WF is always positive. Note 
that for classical CTRWs the limiting phase space distribution 
is uniform and equals l/N 2 both for odd and for even TV. 
Thus, what distinguishes the limiting WF from the classical 
limiting distribution are a few exceptional points in the phase 
space. 

The long time average of the marginal distribution 
|^j(a;;t)| 2 is also recovered from the limiting WF as 
X) K VVj- (x,k) = Xxj, thus for odd N 



(2N-1)/N 2 foTX = j 
{N-l)/N 2 else 



(20) 



t = 1630 
t = 1640 
t = 1650 




100 



40 60 
node x 



100 



FIG. 5: Transition probability (x; t) | 2 as function of x of CTQWs 
on cycles of length (a) N = 101 at times t = 1620, 1640, 1650 and 
(b) N = 100 at times t = 1610, 1620, 1630. The initial node is 
j = 50 in both cases. 

Figure |5] shows the transition probability \ipj(x;t)\ 2 for 
N = 101 and N = 100 at three different times t close to 
the corresponding t r . For N = 101 (N = 100) at t = 1640 
(t = 1620) the CTQW shows the closest revival to the initial 
distribution \i/jj(x; 0)| 2 in that time interval. The correspond- 
ing phase space pictures are given in Fig|5] 

Although the phase space pictures in Fig|6]are vastly differ- 
ent from the ones at the initial time, there are strong peaks in 
the WFs seen as red regions close to the initial point, j = 50. 
Moreover, the patterns in Figs|6ja) and (b) are more regular 
than, say, the ones in Figs|5If) and|4jf). Especially along lines 
perpendicular to the initial node j there are in both cases ar- 
eas with strong positive peaks in the WF, giving rise to strong 
peaks in the marginal distribution J2 K Wj(x, k; t). All this 
means that close to a partial revival the structure is by far more 
regular than at an arbitrary moment in time. 
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V. CONCLUSIONS 



In conclusion, we have defined and calculated the discrete 
WF of a CTQW on a discrete cycle of arbitrary length N with 
PBC. Different from previous attempts, our definition is in- 
dependent of whether N is even or odd. Integrating the WF 
along lines in phase space gives the correct marginal distribu- 
tions. The WF exhibits characteristic patterns. Furthermore, 
we also showed how partial revivals of the WF manifest them- 
selves in phase space. The patterns in phase space allow us to 
draw a much richer picture of the CTQW than the marginal 
distributions can. 



Acknowledgments 



FIG. 6: (Color online) WFs of a CTQW on a cycle of length (a) 
N = 101 at time t = 1640 and (b) TV = 100 at time t = 1620. 
The initial node is j — 50 in both cases. Red regions denote positive 
values, blue regions negative values and white regions values close 
to Oof Wj{x,n;t). 



This work was supported by a grant from the Ministry of 
Science, Research and the Arts of Baden- Wiirttemberg (AZ: 
24-7532.23-11-11/1). Further support from the Deutsche 
Forschungsgemeinschaft (DFG) and the Fonds der Chemis- 
chen Industrie is gratefully acknowledged. 



[1] E. P. Wigner, Phys. Rev. 40, 749 (1932). 

[2] M. Hillery, R. F. O'Connell, M. O. Scully, and E. P. Wigner, 

Phys. Rep. 106, 121 (1984). 
[3] W. P. Schleich, Quantum Optics in Phase Space (Wiley-VCH, 

Berlin, 2001). 

[4] L. Mandel and E. Wolf, Optical Coherence and Quantum Op- 
tics (Cambridge University Press, Cambridge, England, 1995). 

[5] N. C. Kluksdahl, A. M. Kriman, D. K. Ferry, and C. Ringhofer, 
Phys. Rev. B 39, 7720 (1989). 

[6] F. A. Buot, Phys. Rep. 234, 73 (1993). 

[7] P. Bordone, M. Pascoli, R. Brunetti, A. Bertoni, C. Jacoboni, 

and A. Abramo, Phys. Rev. B 59, 3060 (1999). 
[8] W. K. Wootters, Ann. Phys. 176, 1 (1987). 
[9] O. Cohendet, P. Combe, M. Sirugue, and M. Sirugue-Collin, J. 
Phys. A 21, 2875 (1988). 
[10] U. Leonhardt, Phys. Rev. Lett. 74, 4101 (1995). 
[11] A. Takami, T. Hashimoto, M. Horibe, and A. Hayashi, Phys. 

Rev. A 64, 032114 (2001). 
[12] C. Miquel, J. P. Paz, and M. Saraceno, Phys. Rev. A 65, 062309 
(2002). 

[13] J. Kempe, Contemporary Physics 44, 307 (2003). 



[14] R. Iwanow, D. A. May-Arrioja, D. N. Christodoulides, G. I. 
Stegeman, Y. Min, and W. Sohler, Phys. Rev. Lett. 95, 053902 
(2005). 

[15] D. K. Wojcik and J. R. Dorfman, Phys. Rev. E 66, 036110 

(2002) . 

[16] D. K. Wojcik and J. R. Dorfman, Phys. Rev. Lett. 90, 230602 

(2003) . 

[17] O. Miilken and A. Blumen, Phys. Rev. E 71, 036128 (2005). 
[18] E. Farhi and S. Gutmann, Phys. Rev. A 58, 915 (1998). 
[19] O. Miilken and A. Blumen, Phys. Rev. E 71, 016101 (2005). 
[20] U. Leonhardt, Phys. Rev. A 53, 2998 (1996). 
[21] J. M. Ziman, Principles of the Theory of Solids (Cambridge 

University Press, Cambridge, England, 1972). 
[22] O. M. Friesch, I. Marzoli, and W. P. Schleich, New J. Phys. 2, 

4.1 (2000). 

[23] M. Abramowitz and I. A. Stegun, eds., Handbook of Mathemat- 
ical Functions (Dover, New York, 1972). 

[24] D. Aharonov, A. Ambainis, J. Kempe, and U. Vazirani, in 
Proceedings of ACM Symposium on Theory of Computation 
(STOC'01) (ACM Press, New York, 2001), p. 50. 



